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Courses on undergraduate quantum mechanics usually focus on solutions of the Schrodinger equa- 
tion for several simple one-dimensional examples. When the notion of a Hilbert space is introduced 
only academic examples are used, such as the matrix representation of Dirac's raising and lowering 
operators or the angular momentum operators. We introduce some of the same one-dimensional 
examples as matrix diagonalization problems, with a basis that consists of the infinite set of square 
well eigenfunctions. Undergraduate students are well equipped to handle such problems in familiar 
contexts. We pay special attention to the one-dimensional harmonic oscillator. This paper should 
equip students to obtain the low lying bound states of any one-dimensional short range potential. 



PACS numbers: 



I. INTRODUCTION 



A standard undergraduate text in quantum 
mechanics^ devotes most of its attention to analyt- 
ical solutions of the Schrodinger equation. Some texts^ 
coordinate problem solving in quantum mechanics with 
numerical solutions. However, the focus remains the 
solution of the Schrodinger equation. 

A typical introductory course also includes a treatment 
of the more formal part of quantum mechanics, with top- 
ics such as vector spaces and the matrix representation. 
This part is usually more abstract. Examples that are 
often used in this part are two state systems and matrix 
representations of angular momentum. Alternatively, an 
actual (say, 3x3) matrix is given, and students are asked 
to use matrix algebra to mimic a solution of a quantum 
mechanical problem. These examples often strike stu- 
dents as artificial. The main part of such a course appears 
to students to be the solution of various differential equa- 
tions, and the formal part of the course often appears to 
be superfluous. 

Numerical solutions of the time-independent 
Schrodinger equation have been discussed in vari- 
ous contexts (see, for example, Ref. 0). These solutions 
utilize numerical techniques for solving a differential 
equation. Even in the rare instance that an actual basis 
set is used£ it is used to determine the time evolution for 
the time-dependent Schrodinger equation. Many quan- 
tum mechanical research problems that are amenable 
to solution (for example, the behavior of electrons on 
a small lattice) are solved using matrix mechanics. 6 
The purpose of this note is to bring attention to a 
wide variety of time-independent problems that can be 
tackled through matrix mechanics at an introductory 
level. 

In the following section we set up the infinite square 
well problem. The solution to this problem constitutes 
the (infinite) basis set with which we can tackle any sin- 
gle particle potential problem that can be embedded in 
the infinite square well. In Sec. Ill we illustrate the ma- 
trix solutions for the harmonic oscillator potential. This 
example is chosen because the students will have already 



encountered it using the algebraic method (series solu- 
tion) and Dirac's operator method. The former method 
tends to leave students struggling with the mathematics, 
and the latter method leaves many students awestruck. 
The method proposed here uses matrix algebra, which is 
a more familiar mathematical method. 

We also illustrate wave functions, which can be read- 
ily calculated from the eigenvector solutions. We have 
been unable to find problems of this nature in the liter- 
ature. One reference illustrates a solution to the three- 
dimensional harmonic oscillator bounded by two impen- 
etrable walls, but does so following the standard route of 
differential equations, though in cylindrical coordinates J 

The matrix methodology discussed here is at the heart 
of many applications in physics. For example, vibrations 
in solids^ are treated using springs that obey harmonic 
oscillator potentials. Bose condensation in alkali gases^ 
utilizes a harmonic trap, and the conventional theory 
of superconductivity uses electron-ion interactions with 
harmonic oscillators^ Calculations of a single electron 
interacting with vibrating ions (that is, a polaron) that 
use matrix diagonalization can be found in Ref. Illl 

The convergence criterion for the harmonic oscillator 
potential is related to energy scales, and the number of 
basis states required to accurately represent these solu- 
tions is dictated by the energies of the basis states. In 
Sec. IV we briefly discuss some other familiar interac- 
tions, such as the finite square well and the delta func- 
tion potential. In some of these cases length scales also 
play an important role, which leads to a discussion of the 
convergence criterion connected with capturing the short 
length scales in the potential. 

This paper is meant to be sufficiently self-contained 
and should serve to enhance the undergraduate curricu- 
lum in quantum mechanics, and equip undergraduates to 
solve a variety of useful problems. 



II. THE FORMALISM 

A typical introductory course begins with one- 
dimensional problems, including the free particle, the in- 
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finite square well, and the harmonic oscillator potential. 
The free particle has certain pathologies associated with 
the normalization of the stationary states. For this rea- 
son we prefer to discuss the infinite square well, in which 
a particle is confined to the space enclosed by the walls 
of the well. The harmonic oscillator is also amenable to 
analytical solution. Unlike the infinite square well the 
solutions (Hermite polynomials enveloped by Gaussians) 
are unfamiliar to the novice student, and, as far as a 
typical student is concerned, might as well be numerical. 

For these reasons we focus on the infinite square well, 
which is defined as 



Vinf(x) 



if < x < a, 
oo otherwise. 



The Hamiltonian is given by 

h 2 d 2 



H 



2m dx 



2 +V inf (x), 



(1) 



(2) 



where m is the mass of the particle. The eigenstates are 
well known: 
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(3) 



otherwise, 



with eigenvalues, 
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The quantum number n — 1, 2, 3, . . . takes on a positive 
integer value. Students are generally unaware of pertur- 
bation theory at this stage, but reasons for the super- 
script (0) can be readily explained. 

Although we have presented the well-known solutions 
for the infinite square well problem, the emphasis in the 
preceding paragraphs should be on the formulation and 
presentation of a convenient and familiar set of basis 
functions. We are now prepared to tackle a variety of 
problems, represented by H = Hq + V(x), where Hq is 
as in Eq. ([2]), and V(x) is any potential in the domain 
< x < a. We proceed as in most textbooks, except now 
we have the specific basis set Eq. ([3]) in mind. 

We start with the Schrodinger equation in ket notation: 



(H + v)\i>) = E\- 



(5) 



If we use the general expansion of a wave function in 
terms of a complete set of basis states, 



\ip) = ^2 CmlVVn), 



we obtain 



^ c m (H Q + V)\ip m ) = E } j c m \ip n 



(6) 



(7) 



By taking the inner product of ([7]) with the bra (ip n | , we 
obtain the matrix equation: 



(8) 



where 



S E {Q) + 



dxsm( )V(x)sin( ) (9) 



is the Hamiltonian matrix; in the second term we have 
used the eigenstates of Ho as the basis set. As usual, 6 nm 
is the Kronecker delta function. 

Equation ([5]) along with Eq. © allows us to study a 
variety of interesting problems, albeit in the space de- 
termined by the infinite square well. Note that for nu- 
merical solutions of Eq. (jSJ) we need to truncate the sum 
and hence the number of eigenstates used to describe the 
solutions. We denote the maximum basis set number by 
N. 



III. THE HARMONIC OSCILLATOR 

We imbed the harmonic oscillator potential, Vho = 
mu 2 i 2 /2, in the infinite square well. In the following 
all units of distance will be in terms of the square well 
width a, and all units of energy will be in terms of the 
(4) (unperturbed) ground state, Ef^ . We will use lower case 



letters to denote dimensionless energies. The potential 
Vho can be written in terms of the infinite square well 
length and energy scales as 
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(10) 



so that the dimensionless parameter Klo/e[ determines 
the stiffness of the harmonic oscillator potential. We ex- 
pect that for low energy states (say, the ground state), 
the solution should be identical to that of the harmonic 
potential alone, because the wave function will be suffi- 
ciently restricted to the central region of the harmonic 
oscillator potential so that it will not "feel" the walls of 
the infinite square well. High energy states will not be 
well described by the harmonic oscillator results, because 
they will be primarily governed by the infinite square 
well. 

We first use Eq. (J9j) with the potential given by 
Eq. (fini). The result is 
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where 



fjn 



(-1) 



+ 1 



)((n- 



(n — m) 2 (n + m)'* 



(12) 



Note that the g nm remain of order unity close to the di- 
agonal, but for large n the diagonal elements grow as n 2 , 
so the off-diagonal elements become negligible in com- 
parison. 

Equation (TTTj) can be evaluated up to some cutoff for 
a given Uuj/E^; these form the elements of a matrix 
to be fed into an eigenvalue/eigenvector solver (see for 
example, Numerical Recipes^) or in software packages 
such as Matlab, Mathematica, or Maple. 

A typical result is shown in Fig.[T] The low lying states 
are well described by the harmonic oscillator eigenvalues. 
The harmonic oscillator potential becomes truncated by 
the infinite square well at an energy given by the value 
of the potential at x = a (see Eq. [T0|) : i? cross /£^°' ) = 
(■K 2 /16)(huj/E[ 0} ) 2 . For this example, Tiuo/E^ = 50, 
so this crossover occurs at E cmss /E^ w 1550, which 
is approximately where the numerical results in Fig. [T] 
begin to deviate from the analytical harmonic oscillator 
result. For higher energies the numerical results begin 
to resemble those of the infinite square well, albeit with 
an added constant, v . Inspection of Eq. ([9]) (see also 
Eq. (TTTj) ') shows that this constant is the average of the 
harmonic potential. 

What happens if the harmonic potential is consider- 
ably less stiff? If, for example, Tiuj/E^ — 1, then none 
of the eigenstates will describe harmonic oscillator states; 
rather they all will be minor perturbations of the infinite 
square well eigenstates. Nonetheless, to describe such 
an oscillator is still possible, simply by choosing a much 
wider infinite square well. 

How well is the wave function reproduced by this trun- 
cated diagonalization technique? We use Eq. ([6]) and ob- 
tain 

. . , . . 12 . mnx , , 

ip{x) = (x\ip) = > c m \ - sin , (13) 

' V a a 

m— 1 

where the c m are the components of the eigenvectors in 
the basis of the infinite square well eigcnfunctions. For 
numerical results we use Eq. (TIB")) with the coefficients 
obtained from the eigenvectors determined by the matrix 
diagonalization. For the harmonic oscillator the results 
are known analytically. For example, in dimensionless 
units the ground state wave function is 



^ho{x) 



2a 2 E [°) 



1/4 



exp 



(0) 



¥ 

2> 



(14) 

In Fig. [2] we plot the numerical results superimposed with 
the analytical ones. The agreement is excellent (beyond 
what can be discerned on a graph with this scale). We 
can also check the components c m (see Eq. (O) di- 
rectly, because they are determined through diagonaliza- 
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FIG. 1: (color online) Eigenvalues obtained by the exact di- 
agonalization of the harmonic oscillator embedded in an in- 
finite square well (shown as an inset). The symbols denote 
the numerically converged result, obtained with a truncation 
N = 400. We also show the known analytical result for a har- 
monic oscillator, -Eho = hui(n~ 1/2) (we use an enumeration 
that starts with one (not zero) consistent with the enumera- 
tion of the infinite square well eigenvalues). The agreement 
is excellent for low- lying eigenvalues. For eigenvalues that are 
comparable to or exceed the crossover energy the analytical 
results for the infinite square well, suitably modified to take 
approximate account of the harmonic oscillator potential, are 
extremely accurate (dotted curve). We used huj/E[ — 50. 



tion, and compare them with the expected Fourier com- 
ponents of the Gaussian wave function. For example, for 
the ground state, these are 

exp [— n 2 E[ ^ /hw] for n odd, 
for n even. 

(15) 

The exponential decay with quantum number n helps 
to explain why so few components are required to obtain 
an accurate representation of the Gaussian ground state 
wave function. When we checked the numerically deter- 
mined coefficients against those given by Eq. (fT5| (not 
shown), we found excellent agreement. 

The important point is that low lying states for any po- 
tential can be obtained to any desired precision by setting 
up the matrix in Eq. ((9|). Analytical evaluations of the 
integrals are not even required^ For the harmonic os- 
cillator this method requires the least amount of mathe- 
matics. More importantly, it serves to illustrate in a very 
concrete and elementary way the matrix formalism that 
is often presented only abstractly in textbooks. 
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FIG. 2: (color online) The ground state of the harmonic os- 
cillator using the analytical result, Eq. (|14p (symbols), along 
with numerical results, obtained from matrix diagonalization, 
with various truncation values. Convergence is achieved for 
N = 20. We used Tilu/E^ = 50. 



IV. RESULTS FOR SOME STANDARD 
ONE-DIMENSIONAL POTENTIALS 

There are a variety of even simpler one-dimensional 
potentials, and this section illustrates how results from 
matrix mechanics can make contact with known analyt- 
ical solutions. 

Consider the potential of the form, 



V (x) 



if < x < b , 

Vq if b < x < b\ , 

if bi < x < a, 

co otherwise. 



(16) 



This potential represents a generic barrier (well) for Vq 
greater (less) than zero. This case can be solved analyt- 
ically although the algebra is tedious. The units are as 
before, and we let po = bo/ a, p\ = b\/a, and p = p\ — po- 
Evaluation of Eq. ([9|) yields 



TT I pv 



(0) 



v [p- 



+ v (l - 8 nrn ) 



sin 2impi — sin 2nnpo - 
2im 

n — m)npi — sin (n — rn)npo 



7r(n — to) 
sin (n + m)irpi — sin (n + m)irpo 
Tr(n + to) 



17) 



As with the harmonic oscillator, Eq. (fTT)) can be evalu- 
ated up to some cutoff, and the resulting matrix problem 
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FIG. 3: (color online) Eigenvalues plotted as a function of 
quantum number n for the potential shown in the inset. Sym- 
bols denote numerically determined eigenvalues (for a matrix 
truncation iV = 400, but a much smaller value yields identi- 
cal results) . The dashed curve illustrates the analytical result 
expected for an infinite square well of width a/2. The short- 
dashed curve shows the same result for a well of width a. 
These analytical results agree with the numerical results for 
low and high eigenvalues, respectively. Agreement for high 
eigenvalues improves significantly when a shift of «o/2 (dot- 
ted curve) is included to account approximately (as far as the 
high energy states are concerned) for the presence of the step 
potential. We used no = 4000; the value of N is not impor- 
tant if it exceeds the maximum quantum number displayed 
by about an additional 10. 



represented by Eq. ([5]) can be substituted into an eigen- 
value/eigenvector solver. We now consider some special 
cases. 



A. Step potential 

One of the simplest potentials is that of a potential 
barrier as depicted in the insert of Fig. [3] The expres- 
sion in Eq. (fT7|) applies, with bo = a/2 and b\ = a. Fig- 
ure [3] shows the eigenvalues as a function of the quantum 
number n. The symbols show the numerically obtained 
eigenvalues, and the dashed curve shows the result ex- 
pected from Eq. (HJ), with width a/2 - hence the extra 
factor of 4. For low-lying eigenvalues the agreement is 
very good. In other words, the particle finds itself in an 
essentially infinite square well of width a/2. However, for 
sufficiently high energy states the particle should act like 
it is in an infinite square well of width a. We have also 
plotted the result expected for an infinite square well of 
width a, given by the short-dashed curve. Some accom- 
modation should be made for the fact that a step in the 
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potential exists at low energies. (Imagine a state with 
(dimensionless) energy more than 100000 on the scale of 
Fig. [3] — the step with height v = 4000 would appear 
to be a small blip.) The average of this potential is just 
the product of the height and the width, so the dotted 
curve corresponds to the dashed curve shifted by vq/2. 
There are a number of states whose energies are close 
to Wo, where the eigenvalues are not simply related to 
those of an infinite square well. Beyond quantum num- 
ber n = 80 the eigenvalues, with a simple accounting for 
the step Vq/2, agree with those of an infinite square well. 
This example also illustrates how the numerical solution 
converges to known analytical results at low and high 
energies. 



B. Various square wells 

Consider a square well with arbitrary width less than 
a, with Vb < 0. This interaction can be solved ana- 
lytically, although the algebra is tedious. For Vq < 
there will exist states that are bound within the inner 
finite well. Insofar as these states don't "feel" the walls 
of the infinite square well at x/a = 0,1, these energies 
and their eigenfunctions should be given by the simple 
analytical solution found in textbooks. The ground state 
is obtained by finding the lowest positive solution to the 
equation,— 




(18) 



where 6 = b± — bo, u 



E-V 



and uq 



Vq 
<<>> 



a v «r a V e; 

(for Vq < and Vq < E < 0). In the dimensionless units 
used in Sec. IIV Al u = p\/e — vq, and similarly for uq. 

For the finite square well problem analytical results are 
provided by the solution to Eq. (Tig)) , 1 ! and numerical re- 
sults are obtained by diagonalizing the N x N matrix 
with the elements given by Eq. (fT7|) . Specifically, we 
use b = a/2 - 6/2, b x = a/2 + 6/2, and V = V 00 a/b. 
For a fixed value of Poo, this choice of potential (see 
Fig. 2]) allows a study of the eigenenergies and eigen- 
states of square wells with width 6 (0 < 6 < a), keeping 
the product bVo fixed. The limit 6 — > describes an 
attractive delta function potential in the middle of the 
infinite square well, which is treated analytically in the 
Appendix. The other limit, 6 — > a, lowers the baseline of 
the infinite square well by an amount Vq, so that all the 



E n /E 1 



(0) 



V oq /e[ 0) (for 



eigenvalues are shifted: 
6 = a only). 

Results for the ground state energy as a function of the 
finite square well width 6 are shown in Fig. [5J for various 
values of the cutoff A^ (see Fig. Sfb)). The results are 
insensitive to the cutoff for 6 ~ a, while for low values of 
6, there is considerable variation. The inset in Fig. 01a) 
clarifies the situation for 6 w a. The result converges very 
rapidly (just the diagonal element is exactly correct for 





FIG. 4: (color online) Ground state energy as a function of 
finite square well width. The finite square well is imbedded 
symmetrically within the infinite square well, as shown in the 
insert. We use Vbo = —5e[°\ and vary Vb = Vooa/b as we vary 
b. Part (a) shows that convergence as a function of truncation 
(see legend in part (b)) is slow for small b/a. The inset shows 
the situation at large b/a, where the results have converged, 
but not to the analytical result for a square well in an infinite 
space, but rather, as expected, to the trivial result expected 
for a slightly deeper infinite square well. In part (b) we expand 
the low b/a portion of (a), and also show analytical results for 
the <5-function potential, shown as symbols (see Appendix). 
Convergence is slow because of the need to capture very small 
length scales. 



6 




FIG. 5: (color online) Ground state energy wave function 
for a 5-function potential, as calculated by diagonalization of 
matrices truncated at N = 5, 20, 50, and 200. The exact 
result (see Appendix) is given by symbols. There is a cusp 
at the center, which requires smaller and smaller wavelength 
components to capture exactly. We used a 5-function with 



strength bVoo = —5aE\ 
studied in Fig. [4] 



(0) 



corresponding to the limiting case 



b = a). Even the 3x3 matrix yields accurate ground state 
energy values all the way down to b ss 0.5a, and as the full 
frame of Fig.[4ja) shows, finite N results "peel off" as the 
width decreases. Note that the analytical result given by 
Eq. (fl"8|) differs from the numerical result near b = a (see 
inset), because in this regime the two approaches solve 
different problems. The analytical result is for a finite 
square well of width b with no background potential, and 
the numerical result includes the infinite square well in 
addition^ Not surprisingly, the latter case reproduces 
the analytical shift noted earlier from Eq/E^ = 1 to 



Eo/E\ 



(0) 



1 



-4 (with Vqo = -5E\ 



(0K 



Figure HJb) clarifies the behavior in the very narrow 
well regime. Results for the delta function potential are 
given with symbols along the 6 = line. The lowest curve 
is given by the solution of a simple equation (see Eq. (|A2[) 
in the Appendix). There is a considerable variation with 
the cutoff N in the low b/a regime. For example, even 
with a cutoff of N = 100, meaning that unperturbed 
states with energy as high as 10000 are included, the re- 
sult at b = is off from the N — > oo result by about 
10%. The reason is that as the finite well becomes nar- 
rower and the wave function becomes more confined, the 
latter contains components with very small (of order b) 
spatial variations. These can only be described by basis 
states with a similar strong variation, and such states are 
to be found only among the very high energy states of the 
infinite square wellJ^ We use Eq. fl~3")) and show in Fig. [5] 



the wave functions obtained numerically for the case of a 
delta function potential as we vary the cutoff of the size 
of the matrix retained. It is clear that for low cutoffs, 
only variations on a large length scale can be captured. 
Even for the highest cutoff used here, N — 200, the nu- 
merical result still misses (though only slightly) the cusp 
at the centre of the figure. The same kind of plot can 
be examined for square well potentials of varying width, 
so that one can see directly that the matrix size required 
for convergence increases with narrower well width b (as 
already seen through the energy in Fig. [3}. 



V. SUMMARY 

The purpose of this paper is to provide some exam- 
ples accessible to undergraduate students learning quan- 
tum mechanics for the first time that utilize matrix 
diagonalization. 16 Our goal is to provide a context within 
which a simple basis set can be used. The most straight- 
forward case in one dimension is the infinite square well. 
We worked through some examples that are generally 
used in an introductory course on quantum mechanics, 1 
with an emphasis on the harmonic oscillator. We showed 
that by embedding a (truncated) harmonic oscillator po- 
tential in an infinite square well, we can obtain extremely 
accurate solutions. 

We also illustrated how to calculate wave functions and 
demonstrated that potentials with short length scales of- 
ten require a large number of basis states for an accurate 
description. 

The empowering aspect of this approach is that stu- 
dents can tackle bound state problems for any potential. 
For example, it is simple to find the low-lying states for a 
double well potential. The potential does not even have 
to be defined analytically. 13 
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APPENDIX A: DELTA FUNCTION POTENTIAL 
WELL 

Consider a delta function potential well with strength 
a = Vob located at position ao, inside an infinite square 
well between x — and x = a, that is, 



V(x) 



-bV S(x - ao) (0 < x < a) 



oo 



otherwise, 
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where < clq < a (the symmetric case is for ao = a/2). 
The length scale b is introduced to make contact with 
the finite square well potential discussed in Sec. III. The 
application of matching criteria yields the solution (we 
confine our discussion to the symmetric case only): 



The more familiar solution is obtained by allowing a — > 
oo, in which case k = A/2, and a solution always exists. 
In this problem, with a finite infinite square well width 
a, a solution exists only if 



tf)(x) = A 



sinh(Kir)/ sinh(«a/2) 



(x < a/2) 



sinh(ft(a — x)) / sinh(«;a/2) (x > a/2). 

(Al) 

The bound state energy E is given by the solution of 



tanh(«;a/2) = — , 

A 



(A2) 



V 



E 



(0) 



> 



4 a 



(A4) 



/ 2/71 2 771 

where k = J — —^-E and A = — ^-Vob. Equation (|A2j) 

must be solved numerically. The normalization constant 
is given by 



.4 : 



1 



1 - 



A \2 A 



(A3) 



A challenging exercise for students is to check how well 
Eq. (|A4|) holds for a finite square well of non-zero width 
b. Results for non-symmetrically placed delta functions 
can also be readily obtained, but are not as concise, and 
yield a condition for existence of a bound state which is 
more stringent than Eq. (|A4j) . 
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analytical integrals. The numerical approach allows one to 
solve for bound states of any one-dimensional short range 
potential. 

The analytical solution given by Eq. (|18p provides the solu- 
tion for a finite square well in an infinite space. Results will 
therefore differ from the numerical results as the width of 
the finite square well approaches the width of the infinite 
square well. The Appendix provides an analytical solution 
for the 5-function potential in an infinite square well, which 
differs from the textbook solution of an infinite square well 
in infinite space. 

These difficulties are well known. For example, bosons in a 
harmonic trap experience a short range boson-boson inter- 
action. If we try to expand the many-particle wave function 
in terms of the eigenstates of the harmonic oscillator poten- 
tial, convergence difficulties arise because of the disparity 
in length scales. Our problem illustrates this difficulty at 
the single particle level. I am indebted to Doerte Blume 
and Pierre-Nicholas Roy for discussions on this point. 
Though often referred to as "Heisenberg matrix mechan- 
ics" it is now clear that Born should be credited with this 
technique. See, for example, Nancy Thorndike Greenspan, 
The End of the Certain World: The Life and Science of 
Max Born (Basic Books, Cambridge, MA, 2005). 



